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Abstract 



The classical Volterra model, equipped with the Faddeev-Takhtadjan Poisson bracket pro- 
vides a lattice version of the Virasoro algebra. The Volterra model being integrable, we can 
express the dynamical variables in terms of the so called separated variables. Taking the con- 
tinuum limit of these formulae, we obtain the Virasoro generators written as determinants of 
infinite matrices, the elements of which are constructed with a set of points lying on an infinite 
genus Riemann surface. The coordinates of these points are separated variables for an infinite 
set of Poisson commuting quantities including Lq. The scaling limit of the eigenvector can also 
be calculated explicitly, so that the associated Schroedinger equation is in fact exactly solvable. 
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1 Introduction. 



The relation between integrable systems and conformal field theory has long been recognized P3EJ- 
Although the emphasis has been put rightfully on Baxter Q operator and therefore on Sklyanin's 
separated variables |Hj , to the best of our knowledge there is no explicit expressions of the Virasoro 
generators in terms of these variables. 

We make here a first step in this direction by considering the classical version of this problem. Our 
strategy will be to start with the Volterra model on the lattice [IJ |H] equipped with the Faddeev- 
Takhtadjan 0|Bj Poisson bracket. Since the Volterra model is integrable, we can rewrite everything 
in terms of separated variables. Now, the Faddeev-Takhtadjan bracket goes directly to the Virasoro 
Poisson bracket in the continuum limit, and therefore by taking this limit in the separated variables 
formulae we will obtain the Virasoro generators expressed in terms of separated variables. 

This leads to the following rather new type of formula for the Virasoro generators: 

u(x) = £ L n e 2in ™ =pl + 2d 2 x log det Q(x) + (L - p 2 )5(x) 

Here pq is the zero mode and Poisson commutes with everything, the term (Lq — Pq)S(x) will 
be explained later and the formula for Lq is given in eq. (|66j) . The infinite matrix Q(x) reads 
(k, m G {1, • • • , oo}): 

a _ W k (x)d x E m (x) - d x W k {x)E m (x) 

k — ^ 

with 

Wk(x) = — l-Mfc ~ 5 E m (x) = 2mir sin m-KX (2) 

^k ^k 

The above formula for u(x) is valid on the interval < x < 1, and should be extended outside this 
interval by periodicity (in particular the d(x) term in a Dirac comb). 

The result of this paper is that if the variables Z k ,fi k , have Poisson bracket 3 

{Zk,Z k /} = 0, {Z k ,fj, k /} = 2(Z k - plZ^ 1 )fi k 5 k k', {^k,Hk>} = (3) 

then u(x) does satisfies the Virasoro Poisson bracket: 

{u(x),u(y)} = 4(u(x) + u(y)) S'(x - y) + 28"' (x - y) (4) 

Morever, the variables Z k ,/i k are separated variables for an infinite set of higher commuting quan- 
tities, including Lq. 

Since the separated variables are also the ones which solve the classical inverse problem, the 
Schroedinger equation with the potential u(x) 

{-d 2 x - u(x))ip(x, A) = A 2 ?P(x, A)) 



3 Notice that if we redefine Ak = sj Z\ — p§ , the Poisson bracket becomes a standard quadratic bracket {Ak,/ik} = 
2Afc/ife. However po will then enter the formula for Q(x) and, in this work, we prefer to keep that formula simple at 
the expense of a slightly more complicated Poisson bracket. 
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is exactly solvable, meaning that we have explicit formulae for both the potential u(x) and a basis 
of solutions ip(x,A). Constructing the linear combination which is quasi periodic (the so called 
Bloch waves) introduces an infinite genus Riemann surface. The coefficients in the expression of 
this curve define a complete set of Poisson commuting Hamiltonians including Lq. The separated 
variables are points on this curve. 

The paper is organized as follows. In the first three sections we recall some known facts about 
the Volterra model on the lattice. In particular we recall the formulae expressing the dynamical 
degrees of freedom in terms of the separated variables. 

In section 5 we compute the continuum (scaling) limit of the spectral curve. The result is ea. (|41j) . 
We then show that the Hamiltonians H m in this formula are in involution. Moreover we show that 
the scaling limit of the dynamical divisor still belongs to that curve, and hence define separated 
variables for these Hamiltonians. 

In sections 6, 7 and 8, we compute the scaling limit of the eigenvector of the Lax matrix at each 
point of the spectral curve. The result is rather simple and is given in ea. (|5U|) . We then show 
that the obtained expression does satisfy a second order Schroedinger equation and we compute its 
potential u{x). Finally, we construct the two quasi periodic solutions of that equation, the Bloch 
waves, and recover in this way exactly the same spectral curve as the one obtained in section 5. 

In section 9, we give conditions under which the determinants of the infinite matrices that appeared 
in the previous sections exist. We then perform a few checks in a certain perturbative scheme. In 
section 10 we prepare the ground for the serious calculations coming next. 

In sections 11 and 12 we prove that the potential u(x) does satisfy the Virasoro Poisson bracket. 
An essential use is made of certain quartic relations, proven very much like the Hirota-Sato bilinear 
identities. These identities should be considered as generalizations for r-functions of the quartic 
relations on Riemann's Theta functions. 



2 The Volterra model. 



In this and the following two sections we recall some well known facts about the Volterra model. 
The Volterra model, as an integrable system, was introduced in [3]. It is a restricted version of the 
Toda lattice. We consider a periodic lattice with N + 1 sites, and on each lattice site we attach a 
dynamical variable on which we impose the Faddeev-Takhtadjan Poisson bracket [Jj: 



CLj} — didj (4 — Clj — aj)(Sij + i — (^i+l) ~~ a j+l$i,j+2 + 0-i+l^j,i+2 

This bracket 4 is interesting because taking the continuum limit as 

1 



(5) 



en 



1 + A 2 u(x), A 



(N + l) 



4 In terms of Toda Hamiltonian structures, it is a linear combination of restrictions of the second and fourth Poisson 
brackets. 
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it becomes the Virasoro Poisson bracket eq.(|lfl. For precisely this reason, and in this perspective, 
the lattice model has been extensively studied both at the classical level E] and at the quantum 
level |9l I1UI [TT1 1131 IT!?] . The present paper is one more contribution this series of works. 

The Lax matrix for the Volterra model is defined by: 



L(m) 







at 






a 2 



V°i-1 



\/ a iV-l 






/Ojv 



/a N 




(6) 



/ 



It is well known that TrL n (/i) are in involution with respect to the Poisson bracket eq.(jSJ. Hence 
we have an integrable system on the lattice whose continuum limit is directly related to conformal 
field theory. 



The spectral curve V is defined as usual: 

T : det(L(/i) - A) = 
Expanding the determinant we see that it is of the form: 

+ -t{\) = 
where t(X) is polynomial of degree N + 1. 



t(A) = A- X X N+1 - A' 1 (j2 a i) xN ~' + 



(7) 



(8) 



(9) 



where 

A = ^aia 2 ■ ■ ■ a N+ i 

Assuming N = 2n even, t(A) is an odd polynomial, t(— A) = — t(X), and has exactly n+1 independent 
coefficients. However, in that case, there is one Casimir function K = t(2): 

{t(2), ai } = 0, Vi 

The dimension of phase space is N = 2n and we have exactly n commuting quantities. The genus 
of the curve r is g = N. 

At each point (A, /z) of the spectral curve, we can attach an eigenvector *(A, fj,) = (V'i(x)), i = 
1, • • • ,N + 1, corresponding to the eigenvalue A of L(fi). Explicitly, the equation (L(p) — A)^ = 
reads 



aiV>2 + M V^+iV'Af+i 



A^i 



y/o-i-lipi-l + \/a~iipi+l 
/J-y/dN+llpl + y/aNi'N = X^ N+ i 



(10) 
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We extend the definition of the coefficients aj by periodicity cii+N+i = a«> and introduce a second 
order difference operator T>: 



i+1 

This operator is a discrete version of a Schroedinger operator with periodic potential. Eas.()lU[) are 
then equivalent to: 

(vqr^j = Xipi, with 4> i+ N+i = n4>i (11) 
Therefore, the eigenvector \I/ is a Bloch wave for the difference operator T> with a Bloch momentum 

In the continuum limit, ea.()ll|) becomes the Schroedinger equation. 

;>2 „.^™^^„/./'„^ a2„i./„\ \ ~, o a2a2 



(-^-«(iMi)=A%), A~2-A^ 



3 The free case. 



Since in the continuum limit a, — ► 1, it is useful to first recall some formulae in the trivial case 
Oi = 1. They will be generalized to the full case in the next section. To introduce the zero mode 
from the start, we consider the slightly more general case a\ = a: 

\/a(ip2 + (jT^n+i) = Xipi 
■\/a('ipi- 1 + ip i+ x) = Xipi 
y/adiipi + lp N ) = Xip N+1 

The solution of the bulk equations is tpi = ax 1 , + f3x % _ where x± are solutions of the equation 



x 2 -zx + l = 0, x±(X) = -(z±\/z 2 -4), 



X 



Imposing the two boundary equations, we get 

- n) a + (x N+1 - n) (3 = 
- fi)x+ a + (x* +1 - fj,)x- P = 

The compatibility of this system yields the spectral curve 

fi + = x+ +1 + x^ +1 = i(A) (12) 

We now impose that the curve passes through the point (A = 2, fj, = /io) where fiQ is related to the 
value of the Casimir function by 

K = no + Mo 1 (13) 

Setting 

\ a V a 
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ea.()12|) gives a = ijjqr[- Hence the constant a is related to the value of the zero mode po by: 

1 



cos N+l 

The components of the eigenvector, properly normalized, are meromorphic functions on the spectral 
curve ea. ()12j) . Choosing the normalization tpN+i = M> they read 

, r\ \ PN+i-i{z) + fiPi(z) A 

A(\n = 5 — ' z = ~r ( 14 ) 

We have introduced the polynomials of degree j — 1: 

P j (z) = x3+ ~ XJ ~ , P J (z) = zi~ 1 + 0(zi~ 3 ) (15) 

X+ — X- 

The first few polynomials are 

P = 0, Pi = l, P 2 = z, P 3 = z 2 -1, Pa = z 3 -2z 
They are essentially the Tchebitchev polynomials of the second kind. 

As we will see, ea. (|14j) is the general form of the meromorphic function ipi(X, fi) even when Oj ^ a. 
In particular, in order to take the continuum limit, the poles of the eigenvector 5 will have to be 
close to the roots of the equation P N+1 (-^=) = 0, that is 

A(° ) = 2 v ^cos ? ^ T , zf = kv k = l,-..,N (16) 

„r„ l,o,^ ^ . — „±*ATXT or^ c^, /TF71 rrUmv ,,( ) 



For these special values X k , we have x± = e * JV + 1 and ea.(|12[) gives n k = (— 1) . The set of 

points {X k °\ fAiP), k = 1, • • • ,N, will be called the free configuration and will play an important role 
below. 

It is simple to take the continuum limit in this free case. We set 

A~2-A 2 A 2 , z~2-A 2 Z 2 , Vi±i = ^{x ± A), A ' 



N+l 

where we have introduced the variable 



Z = ^V+pl (17) 
The eigenvector equation becomes the Schroedinger equation 

-ip"(x) -pfy(x) = A 2 4>(x), x = Aj (18) 
We also have x± = 1 ± iAZ and the equation of the spectral curve reads: 

M + ^ = (l + iAZ) l ' A + (1 - iAZ)-^ A 

5 In fact in this simple case the eigenvector has no poles at finite distance because they are compensated by zeroes 
in the numerator. This degeneracy is lifted as soon as the en are not all equal. 
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In the limit A — > 0, it becomes 

+ = 2cosZ (19) 
Similarly, the eigenvector becomes a Baker-Akhiezer function: 

, . . sin Z(l — x) + a sin Zx , . 

^(*) = " 20 

sm Z 

When [i = e ±tZ this reduces to ip(x) = e ±tZx as it should be. Notice that when /i is kept as a free 
parameter, the above formula gives two independent solutions of ea. (|18j) . but when [i belongs to 
the spectral curve eq.(|19[) one has 

ip{x + 1) = (J>ip(%) 

Eq. (j20[) presents the two Bloch waves as a single function on the hyperelliptic spectral curve eq. (|19j) . 



Another example, important to us, will be the Dirac comb. 

[-d 2 - H 6(x)]ip(x) = A 2 rp(x), 6(x + 1) = 5{x) 
On each interval Xj = j < x < Xj+i = j + 1, one has 

ip(x) = ctje lAx + f3je~ tAx , Xj < x < Xj + % 

The Bloch condition 

ip(x + xj) = /j, J ^(x), < x < 1 

gives ctj = (fie~' lA yao, f3j = (ne lA y (3o. The continuity of ip(x) gives ao = —^z^yPo, while the 
gap equation on the first derivative, —ip'(xj + 0) + ip'(xj — 0) — Hoip(xj) = 0, gives the spectral 
curve 

i sin A . . 

fi + ^ 1 = 2cosA-H —^- (21) 

The Bloch wave itself is ip(x) = fi J ijj vac (x — Xj), Xj < x < Xj + \, where tp vac (x) is given by eq.(|2()|) 
(with po = 0) but A, [i now belonging to the curve ea.(|21|). 



4 Separated variables. 



In this section, we generalize the previous analysis when aj 7^ a and express the dynamical variables 
of the Volterra model in terms of the separated variables. Equivalent formulae where already 
obtained a long time ago in 0J[3]. A quantum version of this construction for the closed Toda chain 
can be found in |19j . 

We have to reconstruct the eigenvectors ^ of L(n). Let us set ^ = (V>j), i = 1, - ■ ■ , N + 1. We 
normalize the last component tpN+i = Notice that due to eq.Q, [i does not vanish for finite A. 
The components V'i are meromorphic functions on the spectral curve and are uniquely characterized 
by their poles and behavior at infinity which we now describe. 
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We will call P + (A = oo, fi = 00) and P (A = 00,// = 0) the two points above A = 00. In the 
neighbourhood of P^ , the local parameter is A -1 and we have by direct expansion of eq.®: 



P + : fi = A~ 1 \ N+1 (l + 0(\- 2 )) (22) 
P~ : fi = A\- N - 1 (l + 0(\~ 2 )) (23) 



At the points P + and P , the eigenvector Vl/(P) behaves as: 

MP) = 1 = \ l (l + Q(\- 2 )), P~P+ (24) 

^i(P) = >^+i a i«2 • • • a,_i A"* (l + 0(A- 2 )) , P~P- (25) 

This is easily deduced by inspection of eas.l|10|). 

From the general results of the classical inverse scattering theory, we expect g + [N + 1) — 1 = 2N 
poles for the eigenvector (see e.g. [HIE])- Prom ea.()24|). we see that we have a fixed pole of order 
N at P + (on the component ipw), and there remains g = N poles at finite distance, the so called 
dynamical poles. But we notice the symmetry property 

^(-A,-/i) = (-l)>i(A,M) 
so that the dynamical poles come in pairs 

AjV+l-fc = — Afc, HN+l~k = —Hk 

and only (Afe, = 1 • • • n, are independent parameters. 

Everything can be expressed in terms of these 2n = N quantities (A^, ///%), k = 1 • • • n. In fact, they 
can be viewed as coordinates on (an open set of) phase space. 

First, the commuting Hamiltonians are easy to reconstruct. Indeed the spectral curve is determined 
by requiring that it passes through the points (A^,/^), k = 1 ■ • • n, and through the point (2,/jlo), 
where is related to the Casimir function as in eq. (|13j) . 

The equation of the curve itself can be written as a determinant 

/A A 3 • • • A^ 1 n + \ 
2 2 3 ••• 2™ m + 
Ai A? ••• Af +1 w + zxr 1 



det 



(26) 



\ A„ A 3 • • • A^ +1 fi n + fi n 1 / 



Expanding over the first row, we obtain a curve of the form eq.(|8"|). and we can read directly the 
Hamiltonians as the coefficients of t(X). They appear as functions of the (Afc, fik) and can be shown 
to Poisson commute (see ^SJElEi] for a proof and for the quantum generalization of this fact). 
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Eas. (|24l25|) and the data of the N dynamical poles also determine the functions ^ uniquely. Being 
meromorphic functions on a hyper elliptic curve, we can write quite generally 



Q®(\)+fjLR®(\) 

nLi(A 2 -AD 



(27) 



where and i?W are polynomials such that 

Q«(_A) = (-i) i Q«(A), = (— l) i+1 i?W(A) 



Above Afc, we have two points on the curve: (\k,fJ>k) and (Afc, Z^ 1 )- We want the poles to be at 
(Afc,/Ufc) only so that the numerator in ea.(|27|l should vanish at the points (A^,/^ 1 )- This gives n 
conditions 

Q (i) (A*) + (A fc ) = 0, fc = 1 ■ ■ ■ n (28) 

To have a pole of order i at P + and a zero of order i at P~ we must choose 

degree = N -i, degree flW =i-\ 

Hence, these two polynomials depend altogether on n + 1 coefficients which are determined by 
imposing the n conditions eq.(|28|) and requiring that the normalizations coefficients are inverse to 
each other at P^ as in eqs. (|24|25j) . 



It is convenient to use the basis of polynomials Pj(\) given by ea, ()15j) . We will write the formulae 
for ipi in the case i odd, the case i even is similar. 

The polynomial Q (A) can be expanded over 

Q«(A): P 2 (\),P A (\),---P N+1 -i 
and the polynomial R^'(X) can be expanded over 

J?W(A): P 1 (X),P 3 (X),---P l (X) 
Solving the linear system eq. ()28jl . the eigenvector can be written as 



(iiPx (A) ••• fiPi(X) -P N+ i-i(\) 
Fi(Ai) ••• Pi(Ai) -/xiPjv+i-^Ai) 



n(A 2 - a 2 ,) 



det 



Pi (A* 



P(A* 



-^ifcPv+l-i(Afc) 



\Pi(A n ) ••• ^ (A n) ~^nPN+l-i{Xn) 
where .fQ are constants independant of A, /i. Defining 

/Pi(Ai) ••• Pi- 2 (Ai) -MiPjv+i-i(Ai) ••• 



-^(A) \ 
-/UiP 2 (Ai) 

-/x fe P 2 (Afc) 

-// n P 2 (A n ) / 



(29) 



e, 



Pi(Afc) ••• Pj_ 2 (A fc ) -fi k P N+ i-i(X k ) 
\Pi(A„) ••• P_ 2 (A n ) -// n Pjv+i-i(A n ) 



-/iiP 2 (Ai) \ 



-^A,P 2 (Afc) 



(30) 



i^(A„)/ 
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we can compute the leading terms in ea. ()29j) when A — ► oo. At the leading term comes from 
nPi{X), while at P^ it comes from P/v+i-«(A). 

iPi ~ (-1)^ A" 1 Ki det G t \\ P + 
fa ~ (-l)^^detGi +2 A"*, P~ 

Imposing that the two coefficients of A* and A _i are inverse to each other, we get 

ff? = d 

4 det 6j det 6; +2 

Comparing with eas. (|24| l25 [) . we finally obtain 

det 6; det i+3 detB^ . det 9 3 

a N = ~j — «A> o-N+i = — FT"^ V-J-U 



det det 0j + 2 det 0^+2 detGi 

Here A -1 is the coefficient of X N+1 in i(A), eq.©, computed from eq.(|26|). 

We impose the Poisson bracket on the variables Xk,^k 

{Afc,A fc '} = 0, {A fc ,/ifc'} = --S k k' (4A fc - X 3 k ) n k , {/j, k , Hk'} = (32) 

One can then check that the Hamiltonians defined by eq. (|26|) are all in involution (this is a general 
result), and that the Oj defined above do satisfy the Faddeev-Takhtadjan Poisson bracket. The fact 
that the expressions for a7v,a7v+i are different from the ones in the bulk is due to the choice of 
normalization of the eigenvector. However, the Poisson bracket of the <ij is periodic. All this can 
be proved using techniques similar to the ones in 



5 Continuum limit of the spectral curve. 

We now take the continuum limit of the spectral curve ea. ()26|) . The result is ea. ()4ij) . We set 

A = yfaz, X k = y/azk, -7= = 2 cos P ° = z 

yja N + 1 

From these, the scaled variables A, Z, are defined like this: 

A Z Z k 

A = 2cos — , z = 2cos — , Zfc = 2cos— (33) 

N + V N + V N + l v ; 



Notice that we have Z = yj A 2 + Pg . In the following, we will refer to the terminology "perturbation 
theory" when the points (Zk,(J*k) are small deviations from the free configuration ea.(|16|). The 
formulae we will write will make sense in this perturbative setting. This however does not exclude 
the possibility to have a finite number of points which are large deviation. We will also be interested 
in the deviation from the zero mode configuration. That is we make the substitution ^fal — ► \fa^fai 
everywhere on the lattice. Alternatively this amounts to using the variable z = X/y/a. 
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Using the basis of polynomials Pj{z) defined in ea. (|15l) instead of the z 3 , we can write the spectral 
curve as (it has the right form and passes through the right points) 



det 



/ + Pn+2{z) 

Ho + fJ-o 1 Pn+2(zo) 
+ P N+2 (zi) 



Pi{z) Pz{z) \ 

P4(Z ) P 2 (Z ) 

Pa{zi) P 2 (zi) 







XHn+fln 1 PN+2(Zn) P±{z n ) ^(^n)/ 

Without changing the determinant, we can subtract to the first column the linear combination of 
the next two columns: 

PN+2(zk) - Pnizk) = 2 cos Z k 

The first column becomes 

/ fi + /i — 2 cos Z\ 

7o 
7i 

V 7n ) 

where we have set 

lk = Hk + Mfe 1 - 2 cos Z k (34) 
Notice that 70 = 0. The reason for this subtraction is that for the free configuration we also 
have 7^ = 0, so that the spectral curve becomes simply fj, + fi^ 1 = 2cosZ as it should be. The 
subtraction gives sense to the spectral curve in perturbation theory. Expanding the determinant 
over the first row, we can write 

71+1 

H + fJ.' 1 - 2 cos Z = ^2 H2jP2j(z) 

3=1 



The H 2 j are given by 
where we have defined 



H 



I Hn+2 \ 

Hn 



\ H 2 / 



N 



H = N~ X V 

/P N+2 (z ) ■■■ P A {z ) P 2 (zq) 

P N +2^l) ••• P 4 (*i) P 2 { Zl ) 

\PN+2{z n ) ••• P±(Zn) P>i{Zn)' 





/7o\ 


, v = 


7i 


) 


W 



We will need to treat separately the first row and column in the matrix N. Let us write it as 

N -- 



where 



A B 
C D 



P N+ 2(Z ), B = (P N (Z ) ■■■Pi(Zo) P2(Z )) 





(Pn- 


w{zi)\ 














\Pm 


-2(z n ) / 






Pi(zi) P 2 (zi) 
Pi{z n ) PliZn). 
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To zero-th order in perturbation theory, we denote N = and similarly A® , B<® ,C<® , D<® . 

To take the continuum limit we have to consider the matrix NN^ -1 . 



Lemma 1 In the continuum limit, we have 

jVjyW" 1 = f shlZk po sinz fc f i i \q( 1 \m TT7 2„2 ) , fc,m = l,---,oo (35) 

V Z k sinpo Z k \ Zl-mW pg-m^ / ^ ±) ™ TT J 

Proof . Since Pjv+2(<z) = zPn+\{z) — Pn{z) and P^ + i(z^) = 0, we have 

so that , s 

v(o)-i - _JL_ f 1 " BD(0K1 ^ 

A + Bi \ -F (A + Bi)!)^)- 1 -F^BD^- 1 ) 

where F is the column vector with components F k = S k: i, k = 1, ■ ■ ■ , n. It follows that 

/V/v(o)-i _ f 1 

-{^-(C + DF) DD^-^(C + DF)®BD^ 

Noticing that 

A + Bi = P N+2 (z ) + P N (z ) 
(C + DF) k = P N+2 (z k ) + P N (z k ) 

we get in the continuum limit 

1 -{C + DF) k S[nZk Po 



A + Bi ' Z k sinpo 

The main trick to proceed is an explicit formula for the inverse of the matrix Z)(°). It is not difficult 
to check that 

( d,0, -V. = ^w^p 2j (4°>) 

With this, we find 

A 11 

(BD^) m = ^ S m 2 ^£^o)P 2j (4 0) ) 

sm2 #FT 4 -A . 2jp . 2jmir mil 



sin po sin m7r N 4- 1 ^— 
bin JV+1 sin JV+1 iv -f i , =i 

and the last integral is easily evaluated with the result 

fRn (0)-h , 2m 2 7r 2 sinp 

Similarly we compute 

mn (0)-h _ sinZ fc (-ir2m 2 7r 2 
( )km ~ Z k Zl-mW 



E. 2jp . 2jm7T 7717T /" . 
sin— sm— >2 / axsmpnxsmrmrx 
-. . + 1 7v + 1 p y 
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Gathering all this we get ea.(|35[). 



We now introduce the important infinite matrix 

1 



Z| — m 2 7r 2 



and the important vector (77) 



\rf) = M 



-1 



k,m = 1, ■ 
1 

W 



, 00 



(36) 



(37) 



With these notations we can compute the inverse of the matrix NN^' 



Lemma 2 

(jViy(Q)-i)-i = / , v (- ir; ' (-ir JYn 1 \v)(x(po)\ \ M -i\ z k ) 

V sinpo \1-(X(P0)\V) J 2mW ^ m W \^ + P®*J J N± / mfc sinZ fe / 

where we have defined the vector (x(Po)|m = •i_ 1 2 — 2 • 



Proof . With the above notations we can write 

NN^- 1 = 



Letting 



we find 



and 



1 

: Po sin 2 

Z k sinpo Z k 



sinZ,^ sjnZ kMki{dirn _ mxM)2{ _ 1) rn m 2 7r 2 



2m 2 7r 2 [V 1-{X(P0)\V)J ) mk smZ k 
PO ^ 1 N \ ("l) m+1 



( X F) m sinpo \1 - (x(po)|»7) J 2m 2 7r 2 % 



Let us return to the formula for the spectral curve. We assume that the conditions explained in 
Section 9 are satisfied so that the infinite sums we will manipulate are convergent. Denote 

00 

T)(Z) = 1~J2 7 2 — 2 2 (39) 
m=l 

and 

These quantities enter the expression of the continuum limit of the spectral curve. 
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Proposition 1 In the continuum limit, the equation of the spectral curve becomes: 



(j, + \i =2 cos Z + 



sin Z 



Z 



-On 



Z 2 — m 2 7r 2 



where the conserved quantities H m can be taken as 

H m = T m + — (x{Po)\r)Vm, H = Y] 

v(po) 4^ 



H„ 



Pq — 7r 2 m 2 



v(Po 



-y 



Pq — TT 2 m 2 



(41) 



(42) 



Proof . We have 

n n 
H + fl- 1 - 2 C OsZ = Y, P N+2-2i(z)H N+2 -2i =^ / P N +2-2i(z){N- 1 V) i 



i=0 



i=0 



we insert 1 = into the above expression: 



A* + n~ 1 - 2 cos Z = Pn+2-21 (z) (N^- 1 ) ik (N^) kj (N~ 1 V) 3 



Hence, we need to compute 

Y,PN + 2-2i(z)(N^- l ) v 



( P N+2 (z) + P N (z) P N+2 (z) + P N (z) BD ( )- 1 v^ p . 
\P N+2 {z ) + Pn(zqY Piv + 2^o) + Pjv(^) Y ' 



(43) 



whose limit N —* oo is easy to take 



sinZ / p 



2(-l) m mV 



Z \sinpo 

We can now take the limit N — > oo in ea.(|43[) 

sinZ / po 



1 



1 



Z 2 — m 2 ir 2 p^ — m 2 Tr 2 



H + H 1 — 2 cos Z 
where 



Z \ sinpo 



m=l 



1 



1 



/7o\ 



H m = (NN^- 1 )- 1 



71 



V 7 n/ 



x- 



V 



Z 2 — m 2 7r 2 p^ — m 2 7r 2 



71 



v7n 



where X 1 is given in ea.(|38j) and we remembered that 70 = 0. Since 1 — (x{Po)\v) = v(Po) the 
equation of the spectral curve finally becomes 



fJL + fi 1 — 2 cos Z 



sin Z 
Z 



1 



1 



Z 2 — m 2 7r 2 j»Q — m 2 ir 2 



1 



r?7m(x(Po)|r) 
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Another useful expression of this result is: 

M + ^ = 2cc* Z +^fe_i^-4a_Jk_) (44 ) 
Z \ Z z — m A ir z r](po) Pq — m z ir z I 

The next proposition performs a few consistency checks. 

Proposition 2 T/ie points [Z = po, fi^ 1 ), and (Z = Z^,^ 1 ), all belong to the curve eg. \44\ )- 

Proof . When Z = p , we find n + = 2cospo> hence the curve passes through the point 
A = , /if) 1 ) as it should be. When Z = Z k , recalling that f](Z k ) = 0, we find 

M + = 2 cos Z fc + — — }^ 2 _ 2 2 M m; — — ( W + -2cosZ,) 



Z fc ^Z, 2 -m 2 vr 2 m 'sinZi 
2 cos Z fc + S1 ^ fc M km M~j Z{ (in + fj,^ 1 - 2 cos Zi) = n k + ^~ 



Hence the curve passes through the points Z k , /x fe 1 . 



We now show that the H m in eq. 1)41(1 all Poisson commute. We need the following result: 
Lemma 3 One has 

{r n ,r m } = o (45) 

{r n ,r/ m } = {r m ,?7 n } (46) 
{Vn,Vm} = (47) 

Proof . Recall the definitions eas. ()37l4Uj) of rj m and T m . 

Vm=M-l\l) k , T m = M-l\j) k , % = -r-^ lk 

bill Zj k 

where M~, is the inverse of the matrix defined in ea . 1)36)1 . The relation ea.l|47|) is obvious because 
the r] m depend only on the Z k . Consider the second relation ea.l)46|): 

{ Vn ,T m } = {M- k l l k ,M-^ l } = M-]{M- k \^}l k 

= - M ^i M n,U M ^7i}M^\l) k = -M-)M-*{M lp Mrip 

where in the last step we used that {M k i p ,ji} = if k' ^ I. The result is obviously symmetric in 
m and n. Finally the first statement, ea. 1)45)1 . is simple. One has 

{r m ,r n } = {m-i%,m-/7 Z } = -m-^ 
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but because of the structure of M, we have {M rs ,ji} = if r ^ I and for r = I the term in the 
square bracket obviously vanishes. This is a special case of a general theorem |17[ I18j . ■ 



We are now ready to prove 

Proposition 3 The quantities Ho, H m , Poisson commute 

{Ho, H n } = 0, {H n , H m } = 



Proof . Using ea.([42|). one has 

{H n , H m } = Tin({C, T m } + C{C, rim}) - Vm({C, T n } + C{C, r? n }) 
where we denoted 

q — - (x(po)|r) — - ^ 



One has 



hence 



vipo) ^ Po ~ n 2 l 2 



{C,T m } + C{C,r lm } = --±-Cj2 



v(po) i Po ~ 
{T m ,7]i} + {r] m ,Ti} 



All this means that (Z^,^) are separated coordinates for the Hamiltonians H m . 



6 Continuum limit of the eigenvector. 



Having found the continuum limit of the spectral curve, we now consider the limit of the eigenvector. 
Again, the continuum limit can be computed, the result being ea. ()5Uj) . 

As seen from ea. ()29j) . the eigenvector can be written as (for i odd) 

VA det Ni 



U(z 2 - z 2 ) y/det @i det Q i+2 



where 



Ni 



I fiPi(z) + P N+l ^{z) fiP^z) 
Pi{z\) + mP N+1 -i(zi) PiOi) 

Pi(z k ) + fJ,jP N+1 ^i(z k ) Pl(Zk) 



-P N+1 -i(z) 
-HiP N+1 -i(zi) 

-fJ-kPN+l-i(zk) 



\Pi(z n )+fi n P N+ 

1—i \Zn n I 1—i \Zfi) 



~P2(Z) \ 
-/UlP 2 (>l) 

~HkP2{zk) 



(48) 
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Compared to ea.(|29|). we have subtracted the i-th column to the first one for the same reason than 
in the previous section. Also we have used the variable z, Zk instead of A, A^. Let us decompose 
the matrix iVj in blocs particularizing the first row and first column: 

N-( Ui Vi 

where 

Ui = fiPi(z) + P N+1 -i(z) 

(Wi) k = VkPi(zk) + PN+l~i(zk) 

(Vi)j = (j,Pj(z)0(i - j) - P N +i-j(z)9(j - i), i,j odd 



To order zero in perturbation, we have — (— l) k Pj(z^ ) = Pjv+i-j(^ ) so that 



N 



(0) _ (Ui Vi 



o e 



(0) 



where 0^ is the matrix ea.(|3U|) evaluated on the free configuration. It is in fact independent of i 
and we will denote it by 8™. The appearance of zero in the lower left corner was the reason for 
the subtraction in ea.(|48[) and makes things better behaved in perturbation. 

The matrix N- ^ being bloc triangular we can compute its inverse: 



N 



(o)-l 



ur 1 -ur x v&<®- 1 



o e^- 1 

so that 



1 % 



(0)-l 



1 



Returning to the formula for tpi, we multiply all the matrices by 0( o ) 1 . The factors det©(°) 1 
cancel between the numerator and denominator. We arrive at 



jj. det (e^ew- 1 - u^Wi ® y.e^- 1 



det (BieW- 1 I det (©^©(o)- 1 



U k z 2 - 4 



We want to take the scaling limit of this expression. Again, the main trick is an explicit formula 
for the inverse of 0(°). It is not difficult to check that 



K hk N + 1 N+l JV k ' 
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Let us compute 0j0(°) 1 . Using the parametrization ea. ()33|) we find (recall that i is assumed to 
be odd) 



(e.e^- 1 )^ 



mir 



4 s m 7v+T 
iV+lsin^ 



i-2 

E. ]Z k . jrmr 
sm IT? 7 sm TT T 
N + 1 N + 1 

j'=l,odd 

E. (iV + l-j)Z fe . jmvr 
sm — sm 



iV+ 1 iV+ 1 

j=i,odd 

Defining Q(x) as the scaling limit of 0j0(°) _1 , we find (there is a factor 1/2 because the sum is 
over j odd only) 



2 "^~ 1 / dy sin Z k y smrmry - fi k / sin2' fc (l - y) shimmy 
z fc {Jo Jx 



Similarly, we define U(x,A,/j,) and W k (x) by 

17< = (N + (s, A, M ), (Wi)fc = + l)W fe (x) 

We find 

l/(, 1 A,/.)-/.^ + ^^. 2 = v ^s 

and 

sinZ fc x sinZ fc (l-x) 
Wfc(z) = —5; + £ 

^k 

Finally, we have (again there is a factor 1/2 because the sum is over j odd only) 

(v i e^' 1 ) m = v m (x,z, f i) 

where 

, . . f /•* , sinyZ Z" 1 sin(l - y)Z 1 
Vm(x, A, fx) = 2mir i J ay fi — — — smrmry — J — smm7ryj (49) 

Putting everything together, we arrive at (up to a factor 6 independent of x ) 
Proposition 4 

V>(x,A,M) = U(x,A,n) - (V(x, A, n^O^ 1 (x)\W (x)) (50) 

where we denoted by (V(x, A, fi)\ the row vector with components V m (x,A,fj,) and by \W(x)) the 
column vector with components Wk(x). 

6 In this factor we will include in particular ^ ^x-z' 1 /z 2 ) wrucn produces the poles at Z 2 = Z\. This is important 
for the analyticity properties of tp(x) but plays little role for the considerations of this paper. 
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It is easy to show that the infinite sums involved in this formula converge under the conditions 
eq.(|f)5jl of Section 9. 



Ea. (|5U|) is the generalization of ea. (|2UI) , Here, the A and p, dependence is entirely contained in the 
function U(x,A,p>) and the vector V(x, A, p). For the moment they are free complex parameters. 



We now want to specialize to A = 0, /xo = e ±ipo . We have 

U(x,k,n)\ Q ± iPo = _° e ±*Poa: ) V m (x,A,fi)\ Q ± iP0 = U(x,A,n)\ Q ± iP0 V^(x,p ) 



Po 



where 



x,po) = —mir 



+ 



mir ± po mir =F Po 



Hence, up to a constant 

(x,po) = e ±ipox 1 - (V^(x,po)\Q 1 (x)|H / (x)) 

These are the primary fields of CFT. Their logarithmic derivatives are the free fields of the Coulomb 
gaz representation. Notice that we have two such fields playing a completely symmetrical role: we 
go from one to the other by changing po — > —po- This circumstance was recognized and used with 
great profit in |14j . The separated variables make this symmetry explicit and built in. 



7 Schroedinger equation. 

Having found a formula for the wave function -0(x, A,/x), the next question is to find the potential 
in the Schroedinger equation that ip{x,A,p) is expected to satisfy. At this point it is simpler to 
forget the lattice model and work directly with ea. (|5U|) . 

Let us denote by \E(x)) the vector with components 

E m {x) = 2to7t sinrmrx 

Calculating explicitly the integrals in eq. ([49|) . we find (' denotes the derivative with respect to x) 
V m {x,A,p) = U ^ A ^) E ' m (f-n X ^^)E m (x) ^ z 2 = A 2 +p 2 (51) 



and similarly for O(x), we find 



W k (x)E ' m {x) - W' k {x )E m (x) 
Zl — m 2/ K 2 



@km(x) = ^ ' ^2 ' ' (52) 

J k 



Notice the important formulae 

(V'(x,A,fj)\ = U(x,A,fj)(E(x)\, 9'(x) = \W(x))(E(x)\ 
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The derivative of the matrix 0(cc) is a rank one projector. The matrix B(x) has a form familiar 
in the theory of integrable systems, and we know that it leads to non linear differential equations. 
Indeed, let us define the vector |3>(x)) by 

= @- l ( x )\W(x)) (53) 

and let K be the diagonal matrix 

K mm ' TTITT 5 mm ' 



Proposition 5 The vector \<&(x)) satisfies the set of coupled non linear second order differential 
equations 

|*") + 2({§\E))'\<S>) + K 2 \<f>) = (54) 

where 



Proof . By derivation, and using the formula for 0', we get 

\w) = e|$) 

\w') = e'|$) + e|$') = (<s>\e)\w) + e|$'> 

= + 91$") + (($|£>' + ($'|£»|W> 

But we also have 

Z 2 Q-OK 2 = \W)(E'\ - \W')(E\ (55) 
where we defined the diagonal matrix Z 

%kk' = Zu$kk' 

Applying this identity to |3>), we get 

Z 2 G|$) - GK 2 \<S>) = (E'\$)\W) - (E\$)\W') 

But 

Z 2 G\<£) = Z 2 \W) = -\W") 

and therefore 

\W") = (E\<5>)\W) - (E'\<f>)\W) - GK 2 \$) 
Comparing these two expressions of |W"), we find 

6|$") + ({$\E)' + (&\E))\W) = -{E'\$)\W) - ®K 2 \$) 

Multiplying by 6 _1 yields ea. (|5^|) . ■ 

We are now ready to find the Schroedinger equation satisfied by ip. 
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Proposition 6 The function ip{x,A,n) defined by ea. h5U\) satisfies the linear second order differ- 
ential equation 

-t/j"(x,A,fj,)-\p 2 + 2((E\^))'} z/>(z,A,/x) =AV(x,A,/i) (56) 

Proof . We have 

^(x,A lM ) = C/(x,A, M ) - (^(x.A^Jie- 1 !^)) = C/(x,A, M ) - (V(x, A,/x)|$(x)) 
Using the formula for {V(x, A, //)|, we get 

VK^A,/,) = [/ (l - <£'|^-L_|$)) + £/<(i5|_L_|*) 
Next, remembering that U"(x,A,/jl) = — Z 2 U(x, A, //), = — iT 2 |£?(a;)), we obtain 

^(z,A } /i) = -17 (<£|$) + (^|^-L_|$')^ + 17' fl + <^l z2 ^ 2 l^)) 

and 

A, M ) = -C/^ + ^| $ / ) + ((£; | $) y + (jE ;/|_l_|^ 
+[/' ^-(^|$> + (£|^I_|^ 
Using now the equation for |$"), we get ea. (|56(l . ■ 

The potential T{x) = 2({E\Q))' can also be written directly in terms of 0. In fact, we have 
d^logdete = d x Tr e _1 e' = d x {E\Q~ l W) = d x (E\$) 

hence 

T(x) = 2d x (E\§) = 2^1ogdet9(x) 
The Schroedinger equation therefore also reads 

ip"(x, A, n) + [pi + 2 <9^ log det 9] ip(x, A, /i) = -A 2 ip(x, A, /i) (57) 

In this formula both the potential and the function ip{x, A, /x) are known. The potential therefore 
belongs to the class of exactly solvable potentials. It is strongly reminiscent of the formula for 
finite zones potentials |20 [ I21 | 122 } 123] . It can probably also be obtained by an infinite sequence of 
Darboux transformations |24j . 

The parameter \x which enters the function U(x,A,fj,) and the vector {V(x, A, was, up to now, 
a free parameter. Eq. ()50|) therefore provides two linearly independent solutions of eq. ([56j) . We now 
introduce the spectral curve by imposing the quasiperiodicity of ip(x,A,fj,). 
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8 Bloch Waves and Spectral Curve. 



So far, ip(x,A,fj,) was defined on the interval [0, 1]. We extends its definition by imposing 

ip(x + 1, A, fi) = (iif)(x, A, fi) 
This extension is continuous as we now show. 

Proposition 7 

^(l,A,/i)-^(0,A,/i)=0 

Proof . This follows immediately from 

W k (l) = MfcV fc (0), 9 fem (l) = l j.J- 1 Q km (0)(-l) m U(l, A,/i) = MC/(0,A,/x) 



It is worth computing explicitly "0(0, A, /x). In terms of the matrix M introduced in ea.l|36|). we 
have 

e ta (0) = ^(0)M ta ^(0), U(Q,A,fi) = ^, V m (0,A,/j,) = U(0,A,ii) 1 2 2 <(0) 

2^ ^ — 77?. 7T 

It follows that 

where r\ m and are defined in eqs. 1(371 133J1 . Notice that when Z 2 = Z%, we have tp(0,A,fj,) = 0, 
by definition of rj(Z). 

We now turn to the derivative of if/(x, A, fj,). 
Lemma 4 One has 

1/(1, A, /i) - rt''(0, A, m) = /i f (A, m) 
f (A. „) = „ + „- - 2COS Z - ^ ± K ' m '-f - ™^ (59 ) 

z Z 2 — m 2 7T 2 

m=l 



Proof . We have 

sinZ . , . sinZ 

U(l,A,n) = ^—, U(0,A,fj) = — 

and 

A, fi) = fj, cos Z — 1, U'(0, A,fj,)=fj, — cos Z 
Using E k (l) = E k (0) = 0, we get 

^(1, A, M ) = -U(l, A, (J,)(E'(l)\ z2 ]_ R2 + U>(1, A, n) 
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^(0, A, fj) = -17(0, A, M )(£'(0)|-— Lp|$'(0)> + C/'(0, A, M ) 



From this the result follows. 



At this point it is tempting to identify the spectral curve as T(A, fj,) = 0. However, this cannot be 
correct because the point A = 0, \x = does not belong to it. 



We have to change the Schroedinger equation. The only possible modification is at the edges. We 
consider therefore the equation 

i//'(x, A, n) + \vl + 2(EQ)' + H 6(xj\ 4>(x, A, p) = -A 2 ^(x, A, fx) (60) 

The bulk formula for t/j(x,A,n) does not change. The continuity of tp{x,A,fj) at x = 1 still holds, 
but the derivative now has a discontinuity 

r i+ 

J dx jp" + H if>(l) = 

Using 

/ s sin Z 
^\,A^)=nr){Z) — 

the Bloch condition becomes 

A, h) - ^'(0, A, fj) - H fi T){Z) = 

that is 

H + fi- 1 = 2 cos Z + ^ \J2 z2 T _ m m2 ^ ~ H Q -q{Z)\ 

where we have set 



r m = 24(0)34(0) - <(i)$L(i) (6i) 



We now determine the coefficient H$ by requiring that the curve passes through the points po, /^j 
We find 

h = 1 r "t 

Hence the curve takes the form 

_i sinZ ( r m v(z) r m \ 

fJ, + fl = 2cOSZ + — — > =-j 7y~\~2 2~~2 



±1 



In order to compare with ea. (|44j) . we must compute T m . We have 

^•(0) = ^^^, W fc (l) = ^^, Wtt0) = l-n k coiiZ k , W' k (l) = cosZ k -n k 
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and 

9 fcm (0) = W ^ )E '^ = W k (0)M km E' m (0), 6 fem (l) = W ^ l)E '^ = W k (l)M km E' m (l) 
Zj; — TT z m z Z£ — TT z m z 

where M is the matrix introduced in eq.®. Since |$'(0)) = 9" 1 (0)|W(0)) and |$'(1)) 
e- 1 (l)|T^ r '(l)>, one has 



^(0)C(0) = ^^( / x- 1 -cosZ fc )-^, E' m {l)& m (l) = J2M-l(co S Z k > Z/ ' 

sin Z//; 

A; A; 



MfcJ 



sinZ fc ' rov ' mv ' V mfc ™sinZ fc 

hence 

r™ = fL (ciW.. (o) - fL n W- m = V w-J 

sin Zfc 



r m = £4(o)*' m (o) - ^(i)C(i) = E M -fc ^flr 7fe 



where we recall that 7^ = ^ + 1 — 2 cos Z/% . This is exactly eq. (|40[) and shows that we have 
recovered precisely the spectral curve ea . (144(1 . 

Finally, let us close this section by proving the 

Proposition 8 The function ip(x,A,fi) has no pole at the point Z = Z k ,y^ 1 

Proof . Here we restaure the factor „ , 1 a , 2 . . 

U(x,A,n)\ z=Zk!ll -i = n^Wkix),^ V k f(x,A,fi)\ z=Zk ^-i = n k l &w{x) 

t[)(x,A,fj,)\ z=Zk ^-i ~ ^ Z~zyz2j ^k 1 W k (x) - Q kk '(x)e^ 1 kll (x)W k// (x) = regular 



This shows that the same property for the eigenvector on the lattice has been preserved when 
taking the continuum limit. 



9 Perturbation theory. 



In the previous sections, we have manipulated determinants of infinite matrices quite freely. It is 
necessary now to investigate the conditions for the existence of the determinant det 0(x). We recall 
the free configuration ea. (|16|) . In the scaled variables it reads 

4 0) = ^ /4 0) = (-D fc 

By construction, when (Z k ,fi k ) = (Z^\ [J< k ), we have Q(x) = Id, so that detO(a;) = 1. Clearly 
for det0(x) to exist, we have to assume (Z k ,fj, k ) — ► (Z}^\ fJP^) when k — > 00. Hence we set 

Z k = kir + 6Z k , fx k = (-l) k (l + 5fi k ) (62) 
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It is not difficult to see that to leading order in (5Z k , S/Jj k ), we have 

Wk (x) = - — (5Z k cos kirx — 5fi k sin kirx) 

K7T 

and this implies 

2 ?Ti 7T ( 

®km(x) = —, — o n — ^Ti 8Z k (mcos kirxcosrmrx + k sin kirx sin rmrx) (63) 

k[Zi — ir 2 m 2 ) I 

+8fi k (—m sin kirx cos rmrx + k cos kirx sin rmrx) \ 

Notice that when m = k this formula gives that to leading order @k,k(x) = 1, as it should be. 

A first consequence of these formulae is that if 5Z k = 0, 8ji k = 0, beyond a certain index k = k max , 
then for k > k max we have Wk(x) = 0, ®k,k(x) = 1> and &k,m{x) = 0, Vm ^ k. As a result only 
the first block of size k max x k max of the matrix Q(x) plays a role and all the constructions of the 
previous sections reduce to finite size matrices and vectors. 

If however we want to retain an infinite number of modes in order to keep the field theoretical 
character of the model, one has to say something about the rate at which bZ^ and 5^ tend to zero 
when k — > do. Disregarding a finite number of possibly large bZ^^bii^ which play no role in these 
convergence questions, we may assume that @(x) is given by eq. (|63|) . As we have seen, it is of the 
form 

e(x) = id + e(x) 

where Q(x) is small. In fact, bounding the trigonometric functions by 1, we have 

— 77J 

|O fe , m (x)| < c— ?(\5Z k \ + \8[j, k \), m^k 

k\k — m\ 

Since \k — m\ > 1 when k ^ m, we may write as well 

\Qk,m(x)\ <c—(\5Z k \ + \5fj, k \), m^k 

It is not difficult to see that this formula is valid also for m = k (we have to adapt the constant c). 
It follows that 

logdetG(x) = Trlog(l + Q(x)) 

00 i 00 n ( \ n 

< E- T ri (*)r<E^ E(i ffz *i + M ( 64 ) 

n=l n=l \ k / 

Hence a sufficient condition for the existence of the determinant is that the series M5Zfc| + |5/ifc|^ 
converges. This is achieved if 



\5Z k \ + \6ii k \ < e > (65) 

One can then adjust the constant d such that the series in ea. ()64j) converges. The condition ea. (|65() 
ensures that logdetG(x) exists. To build the potential u(x) however, this function has to be twice 
differentiable in x and this may require stronger conditions on e. 
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Now that we have found an expression for the potential u(x) in terms of a countable set of variables 
Zkil^ki we would like to check the Virasoro Poisson bracket directly. Recall that 



„2inx 



u(x)=p 2 + T(x)+H 5(x), T{x) = 2d x {E\$>) = Y J T n e 2 

n 

Notice first that T(x) has no Fourier component To: 

T = / dxT(x) = 2 f dx((E\$))' = 2( (E(1)\§(1)) - <£(0)|$(0)> ) = 
Jo Jo 

where we used that E m (Q) = E m (l) = 0. The Fourier expansion of the potential u(x) reads 

u(x) = £ L ne 2mnx =P 2 + + H )e 2l ™ x 

n n 

we must therefore identify 

L =p + H =p + — — > — — 66) 

and 

L n = T n + H , =>- T n = L n - L + pi, n^O 
If u(x) has Poisson bracket eq.@, the algebra of the L n reads 

{L n , L m } = 8in(n - m)L n+m - 16iTT 3 n 3 5 n+mfi 



{T n ,T m } = 8iir(n — m)T n + m — SmnT n + 8mmT m — 16i(7r 3 n 3 — isp\n)& ri 



The Poisson algebra for the T n is then closed: 

Jn+m,0 

or, in a form that will be useful later 

{T(x),T(y)} = 25"' (x - y) + 4(2pg + T {x) + T(y))8'(x - y) - 4T'{x)8{y) + 4T'(y)5(x) (67) 

In this section, we consider the situation where all the variables {Z^,^) are close to the free 
configuration as in ea. (|62|) . and we perform a perturbation theory in 5Z k ,5^k- We have seen that 
to lowest order ©(°'(a:) = Id by construction. So, we can write the expansion 

0(x) = Id + G (1) (x) + © (2) (x) + ••• 
\W{x)) = \W w {x)} + \W (2) (x)} + --- 

where we have taken into account that |W(°'(a;)) = 0. It follows that 

|$(a?)) = |$( 1 )(x)) + |$( 2 )(x)) + --- 

where 

l$W) = |cd( 2 )) = \w&) _ e«|w«), • • • 

To lowest order, we find easily 

rW(z) = 2{E\¥^)' = 4^ kir (5 Z k cos 2knx - Sfj, k sin 2kirx) 
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This shows in particular that 5Z k and 5/j, k are just the Fourier components of the potential in this 
first approximation. We see here clearly that for T^\x) to exist as a function (and not just as a 
distribution), we need e > 1 in ea.(|65|). The Poisson bracket eq.© becomes to leading order 

{6Z k ,6ii k } = 2hr(l ^ 



fc 2 vr 2 , 

To define modes independent of the zero mode po we introduce 

«fc = a k (5fj, k - iSZ k ), a\ = a k (5fi k + iSZ k ) 
where the coefficients a k , a k satisfy 7 

1 k 2 TT 2 

47TPq — k z ir z 

With this choice, one has 

{a k ,a[} = ik 

we can then rewrite 

T^\x) = 2i Y kit ( ^e 2lk ™ - i e ~ 2lk A 
k W «* J 

It is now straightforward to compute the Poisson bracket 

{T^(x),T^(y)} = 25"'(x-y)+8p 2 5'(x-y) 

This is the correct result for the Virasoro Poisson bracket in this approximation. Notice that the 
term 5"'{x — y) is exact already at this level. Higher order terms cannot contribute to it. 



Next, we look at the conserved quantities. The leading terms in the expansions of r\ m and H m are 
easy to find: 

rim ~ 2mvr 8Z m , T m ~ 2m 2 vr 2 ((5 / u 2 n + SZ^) 
To see it, consider the defining relations of T\ m 

V v2 Vm 2 2 =i,vfc 

^ Z 2 - ir 2 m 2 

m K 

When Z k is given by ea. (|62|l . the dominant term in the above sum is m = k. The equation becomes 
r] k /(2k-K5Z k ) = 1. The same argument starting with the equation of the spectral curve, eq. 1)41(1 . 
taken at the point Z k , /j> k which belongs to it, yields the formula for H m . Remark that the condition 
eq . 1)65)1 ensures that the sums in the definition of the function r](Z), eq. 1)39)1 . or in the definition of 
the spectral curve, eq. 1)41)1 . are convergent. 

Written with the oscillators a m ,aln, we find 

H m = 8tt(pI - ■K 2 m 2 )al n a m , H = 8ir a^ m a m 

m 

7 It is known |14| that the poles at p§ = k 2 ir 2 are classical remnants of the zeroes of the Kac determinant. 
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It is clear that the H m are in involution at this order. As we see, in first approximation, the 
dynamical system reduces to a set of decoupled harmonic oscillators. The generator Lq is given by 

Lq = pi + 8tt ^2 a^fim 

m 

It is easy to verify that 

{L ,T^(x)} = -4d x T (1 \x) 



These perturbative arguments are good indications that u(x) indeed satisfies the Virasoro Poisson 
bracket. Clearly we will not go very far in perturbation and we now look for a more formal proof 
of this fact. For that purpose, we need some preparation. 



10 Some identities. 



Before computing Poisson brackets to check the Virasoro algebra, we collect a number of useful 
identities. We start with a formula for the inverse matrix G~ 1 (x). It has the same form as @(x). 

Proposition 9 Let us define 

(F\ = (EIQ- 1 (68) 



Then, we can write 



Proof . Multiplying eq.(|55|l on both sides by G , we get 
and so 

i (@- 1 \w')) m ((E\&- 1 )k - (e- 1 \w)) m ((E'\e- 1 ) k 

mk 7r 2 m 2 - Z\ 

But 

e -1 ^') = |$'> + (£|$>|$>, {E'\Q- X = (F'\ + (E\$)(F\ 
Plugging into the above formula, we obtain eq. (|69|) . 



Proposition 10 The vector (F\ satisfies a set of differential equations. 

(F"\ +2((E\$))'(F\ + (F\Z 2 = (70) 

Proof . The proof is the same as for |$). 
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From this we easily deduce 

(e-iy = -\$)(F\ 

which can also be proved using the similar property of G. Let us define 

&( \ _ E m {x)$ m {x) , , _ ^ E' m (x)$ m (x) 

Zr — m tt Zi — m z ir z 

m K m K 

{ ~'k\ x ) = / T7P2 2 2\2> u k\x) = 7 TTPl o 2^2 

^ [Zr — m z 7r z ) z *— ' [Zr — m z ir z Y 

m k ' m y k > 

Proposition 11 We have the identity 

(1 - B k )W k + A k W' k = (71) 

Proof . This is just a rewriting of \W) = @\&) using eq. (|52|) for O(x). ■ 

Proposition 12 The following two identities hold 

(l-B k + A' k )F k -A k F' k = (72) 
(B' k + Z 2 k A k )F k + (l-B k )F k = (73) 



Proof . The first identity is a rewriting of (F\ = (i?|0 -1 using eq. (|69ft for _1 (x). The second 
identity is just a rewriting of (F'\ = - (E\$)(F\. ■ 

The above two identities form a linear system for F k and F^. Its compatibility implies the following: 



Proposition 13 

(1 - B k ) 2 + A k B' k - A' k B k + A' k + Z\A\ = (74) 

Let us define 

F k = WUl-B k ) + W k 'A k (75) 
An important consequence of eas. ()72l73|) is 

Proposition 14 The functions F k {x) and F k {x) are proportional. 

F k {x) = --^-F h {x) (76) 
fJ-klk 

The proportionality coefficient is written in this specific way for later convenience. The quantity 7^ 
is the one defined in eg. \34\j- 
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Proof , Let us compute the Wronskian 
Wr(F h , F k ) = F k (w' k \l - B k ) - W' k B' k + Wj!'A k + W^) - F' k (w' k {\ - B k ) + W^A k ) 

= -k 2 k W k ((l -B k + A' k )F k - A k F' k ) - W' k ([B' k + Z 2 k A k )F k + (1 - B k )F' k ) = 



Proposition 15 The following quantities are constants independent of x 

Vm = <5> m E' m - & m E m - <E$> <f> m E m - 7t2{ 2 1 _ m 2 ) (^K ~ ®'m%){ E 'm E V ~ E ™ E ' P ) ( 77 ) 

The quantities rj m defined in eg. \77\) are in fact the same as the ones introduced in eg. 

Proof . To prove the first statement, just take the derivative with respect to x and use eq. (|54j) . To 
prove the second statement, use the 7/ m defined in ea. ((77|) to rewrite ea. ()74j) as 

E (7 2 —2 2 x =l. Vfc ( 78 ) 
^— ' (Zf — m z TT z ) 

which is exactly the same as ea. (|37j) , ■ 

A straightforward consequence is the "trace" formula that will be useful later: 

(E&) - + {E$) 2 = - J2 (79) 



■m 



We now compute the coefficients ( k appearing in ea. (|75|l 

Proposition 16 The coefficients ( k in eg. j76\ ) are determined by the set of equat 



ions 



E 7 2 ^ 2 2 = i> Vm ( 80 ) 

~* Z k — m z TT z 

These equations are dual to eq. &78\) . 
Proof . Start with 

(f\w) = (E\e~ l \w) = (e\<s>) 

Using eq. (fTTl 1751761 ) , we have 

F k W k = ^-(-W? + WZW k )A k = C k A k (81) 

Hklk 

(F\W) = E CkA k = E ( E z z - k mw ) EA = {m = E EA 

k m \ k k J m 

Since this has to hold for all x, the only possibility is eq.(|80|). ■ 

We collect below a few more identities of the type of eq.([81j) that will be important later 



hence 
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Proposition 17 



F k W k = ( k A k 

F k W' k = -C k (l-B k ) 

F' k W k = ( k (l-B k + A' k ) 

F' k W' k = Ck(B' k + Z 2 k A k ) 



Next, we relate the r] m and ( k 

Proposition 18 The following relation holds: 



^ (z 2 - m^y a 

Proof . We start from 

rn 

When k = k! this gives 

W k F' k D k - W k F k {D' k -A k + Z 2 k C k ) - W k F' k C k + W k F k (C' k - D k ) = I 
or using eas. ()82l85|) 

C k {2D k (l - B k ) + D k A' k - D' k A k + A\ - 2Z 2 k C k A k - C k B' k - C' k {l - B k )} = 1 

Expanding this formula using (p ^ m) 

1 1 1 

{Z 2 -Ti 2 m 2 ) 2 {Z 2 -vrV) ~ ~ir 2 (p 2 -m 2 ) (Z 2 - ir 2 m 2 ) 2 

1 f 1 1 



7T 4 (p 2 - m 2 ) 2 1 (Z 2 - 7r 2 m 2 ) (Z 2 - ttV) 
we get the result with rj m represented by ea. ((771) . 

An immediate and important consequence is an expansion of the function r){Z) near Z 

71 (Z) = (Z 2 -Z 2 )±- + ... 

Returning to the formula for ip(x, A,/i), we can write it as 

ip(x,A,n) = - ^ — w(x,Z) - - '' * w*(x,Z) 
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where we define 

Mx,Z) = e^^-iE'l^^m + iEl^^m) (87) 

w*(x,Z) = e- iZ *(l-(E>\z^\*)-(E\z^\*] 

The functions w(x,Z) and w*(x, Z) are defined on the Riemann sphere with a puncture at do. One 
can compute the Wronskian 

w'w* — w*'w = 2iZr](Z) 

This Wronskian vanishes precisely when Z 2 = Z?. Hence, when Z = Z k , w(x, Z k ) becomes 
proportional to w*(x, Z k ). Indeed we have 

w(x,Z k ) = e iZhX {l-B k + iZ k A k ) 
w*(x,Z k ) = e- iZ " x (l - B k - iZ k A k ) 

Then eq,(|71|) can be rewritten as 

a {+) 

w(x,Z k ) = * w *( x ,Z k ) (89) 
a) 



where 

a k ±] = 1 - fi k e ±iZk , 4 +) 4~ } = Vklk (90) 

The function F k (x) introduced in eg . (|75j) is solution of eq.(|7()jl. It is not difficult to see that the 
other solution of this equation is 

G k = A k W' k + 2Z k (W k D k - W' k C k ) + xF k 



In terms of w(x, Z), w*(x, Z) defined in eas. (|87l88j) . we have 

% = ^{a { - ) w\ Zk + a k +) w*\ Zk ) = OL k ~ ] w\z k = a { k +) w*\ Zk (91) 
G k = —(a^dzwlz.-a^dzw*^) (92) 
This will play an important role below. 



11 Virasoro algebra. 



We are now ready to compute the Poisson bracket {T(x), T(y)}. The result is precisely the algebra 
of the T n = L n — L , ea.(|67|). 
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Proposition 19 Let T(x) = 2d x (E(x)\&(x)) , and X(x),Y(x) be two arbitrary test functions. 
Then we have 

if dxX(x)T(x), [ dyY{y)T{y)\ = - [ dx(X"'Y - XY'") - 4 f dx(X'Y - XY')(p 2 + T) 
I JO Jo J Jo Jo 

+4 I dxX(x)5(x) [ dyY(y)T'(y) -4 f dxX(x)T'(x)[ dyY(y)5(y) (93) 
Jo Jo Jo Jo 

The proof is rather long and we will split it into several lemmas. Since |<£) = Q~ 1 \W), we have 

|$) 2 } = e 1 - 1 e 2 1 [{6i, e 2 }|$ 1 )|$ 2 ) - {©i, w 2 }|*i> - {w u e 2 }|$ 2 ) + {w h w 2 } 

where the index 1, 2 refers to the customary tensor notation. In this expression, all Poisson brackets 
can be computed explicitly. Using the fact that rows with different indices in G and W Poisson 
commute and using only ea. (|71|) . we arrive at 

{d> m (x), = -2^e m 1 fc (x)G; fe 1 ( y )^ (l - [F k {x)G k {y) - F k {y)G k {x)] (94) 

where 7^ are defined in ea. ()34|) . Multypliying by E m (x)E n {y) and remembering that (E\Q~ l = (F\ 
we get 

{(E(x)mx)),(E(y)My))} = 2j2^^(l-$) x (A k {y)B k (x) - A k (x)B k (y)) (95) 

k k1 k H k \ k / 

where 

A k {x) = Fi(x), B k {x) = F k {x)G k {x) 

Using eas. ()87l88|) . we can write 

A k {x) = /x fc 7 fc w(x,Z)w*(x,Z)\ z =z k 

B k {x) = f^{w*{x,Z)d z w{x,Z)-w(x,Z)d z w*(x,Z))\ z=Zk 

The strategy to evaluate the right hand side of eq. (|95|) is to rewrite it as a sum over the residues 
of certain poles of a function on the Riemann Z-sphere. This sum can then be transformed as a 
sum over the residues of the other poles (there will be none in our case) plus a integral over a small 
circle at infinity surrounding an essential singularity. This last integral can then be evaluated using 
the known asymptotics of the function. Let us define 

and introduce the functions 

Oi = w(x, Z)w* (y, Z) — w* (x, Z)w(y, Z) 
= a(Z)w(x, Z)w(y, Z) 

+b(Z)(w(x, Z)w*(y, Z) + w*(x, Z)w(y, Z)) 
+c(Z)w*(x,Z)w*{y,Z) 
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These functions are denned on the Riemann Z-sphere have poles at the points ±rmr and have an 
essential singularity at infinity. Let 

and recall the definition of rj(Z) eq. (|39[) . 



Lemma 5 

£ ReS nk) i 1 ~ i) = ~ 2 S tSS f 1 " I) X O 4 ^^ " M*)B k (yj) (97) 

-J-2, I \ J \ / p, k ! k h' , k \ fe / 



Proof . The factor r] 2 (Z) introduces double poles at ±Z k because f]{Z k ) = 0. However using ea.(|89|). 
we see immediately that fii|±z fc = 0, so that the poles are in fact simple. Remembering ea. (|86|) . 
we have 

E^(i-|) = eJ('-})K+Wa) 

We need to compute dz^\±z k — dz^i\±z k ^2\±z k - Evaluating at gives 

2z f (y^~ oc \ / \ 

c^k = -2-2 U(Z fc )-^ + 26(Z fc )+c(Z fc )-^ U fc (y)i3 fc (x)- A(x)i3 fc (y) 

while using w(—Z k ) = w*(Z k ), dzw\-z h = — <9zii>*|z fc5 we also have 
2% ( , „ ,ut „ „ N , „ N a 

The result follows from the identities 

. at , , „ s , «7T Zj. sin Z k 

a(Z k )^ + 2b(Z k ) + c(Z k )^ = 2i- k 



d z ^t\-z k = -r^ c(-Z k )^t + 2b(-Z k ) + a(-Z k )^ ) (Ak(y)B k (x) - A k (x)B k (y) 



Ik 

c(-Z k )^t + 2b{-Z k ) + a{-Z k )^r = 2i 7 " 



a k , nu 7 n , { 7 \ a k n -Z k smZ k 



Ik 



Next we have to examine the poles at ±m-7r in the expression 



r7 2 (Z) V ^ 2 
We rewrite fii and as 

fii = (w(x) — w*(x))w*(y) — (w (y) — w*(y))w*(x) 

n 2 = — Z — x Uw*(x) - w(x))(e lZ w*(y) - e~ lZ w{y)) + 
41 sin Z l 



+(w*(y) ~ w(y))(e iZ w*(x) - e~ lZ w(x))} 
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Recalling the formula ea. (|871 I5S|) for w(x,Z) and w*(x,Z), we see that when Z = 0, we have 
w = w* so that fii = 0(Z) and ft 2 = 0(Z 2 ). Hence we have no pole at Z = 0. When Z = ±7rm + e 

1 (-1) . (o) * 1 *(-l) . *(0) . 

with 

w ±m( x ) = w ±m\ x ) = =Fvrm$ m (x) 

Because the two leading terms are the same, both w*(x, Z) — w(x, Z) and e lZ w*(x, Z) — e~ lZ w(x, Z) 
are regular. So f^i and f2 2 both behaves like 1/e. Since l/r] 2 (Z) behaves like e 2 , the whole thing is 
in fact regular. 



We come to the conclusion that everything happens at infinity. We want to compute 

dxX(x)T(x),J dyY(y)T(y)\ = 4 j\xX{x) dyY{y) dZ (l - ||) -J—^fi 

(98) 

Let 

w(x,Z) = ir 1/2 (Z)w{x,Z) 
w*(x,Z) = r ) - 1 / 2 (Z)w*(x,Z) 

The wronskian of w(x,Z) and w*(x,Z) is 2iZ and therefore these functions coincide with the 
Baker- Akhiezer functions which are usually introduced in the pseudo-differential appproach to the 
KdV hierarchy (see e.g. |15j). At Z ~ 00, we have 

JZx ( , U(x) J{x)+U} 2 {x) 



_,-/,/ , uix) uj'(x) + uj 2 (x) 
w*{x,Z) = e lZx \ + -±-L + '.^ + 



iZ 2(iZ) 
where we have set 

u(x) = (E&)(x) 

Ea. ()79() is needed to verify this formula. Using theses asymptotic forms, we find 

/ 1 \ 

2iZx 



d x {w 2 (x,Z)) = Yl Mx){iZ) n \e 

\n=— 00 / 

d x {w* 2 (x,Z)) = ( £ (-l) n A n (x)(iZ) n j e~ 2t 

\n=— 00 / 

-1 

d x {w(x,Z)w*(x,Z)) = Yl C 2n(x)(iZ) 



\2n 



where 



X A l-X 



Ax =2, A = -4w(x), A_i = 4w 2 (x), A_ 2 = - jV - 2 / u/ 2 , A_ 3 = ^w 4 + 4w / 

C-2(^) = w"(x) 
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Consider the term proportional to b(Z) in ea.(|98j): 
dZ 



4 



^,b{Z) (l- ||J J dxdy(X(x)Y(y)-X(y)Y(x))d x w 2 (x,Z)d y w* 2 (y,Z) 



Lemma 6 Let us define 



We have 

I-i(z) = - l -5(z), I ( z ) = -±6>(z), L l {z) = -U"(z), / 2 (z) = -if( z ) 

on— 3 

Proof . It is clear that 

d z I p (z) = 2I p+ i(z) 

hence we can determine all the L p (z) recursively. For p > this is done by successively differentiating 
I-i(z) which is easy to calculate 

7 _ l(z ) = /" ^ h{ z) -I e 2iZz = - f *?_^L e xzz = -ly e 2 — = -U{z) 
y ' Jr 2z7r v ' iZ 2 7 r 2ivrsinZ 2^ 2 v ; 

For p < —2 we have to successively integrate I—\{z). For this we need boundary conditions which 
are provided by 

/ ^-b(Z)Z~P = 0, p>2 (99) 

This is because 

[ ^Lu Z )Yc?Z~v = -i / ^^y C P-i Z -p+i 
7 r 2ivr v ; ^ s 2i 7 r 2ivr sin Z ^ s 

( 2 f dZcosZ 1 C 2 ( , 1 ^ 2C 

cot c - - - 2^ 



JOoo S \ \> n=Q «5 



Denoting 



F p (x,y) = {-^) m Mx)A m {y) 



we get 



m+n=p 



•1 r l P= 2 



4 / dx [ dy{X{x)Y(y)-X{y)Y(x))y j F p (x,y)(I p (x-y)+pll p - 2 (x-y)) 
Jo Jo 
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In this expression, we separate the terms with a 5(x — y) function or its derivative which will lead 
to local terms (X&) and the non local terms (iVX&) with are proportional to e(x — y). 

L b = 4 J dxdy(X{x)Y(y) - X(y)Y(x)) (f 2 I 2 + F X I X + (F + p 2 F 2 )I + (F_j + p§Fi)I_i 

/oo 
dxdy(X{x)Y(y) - X(y)Y(x)) ^(F_ 2 _ n + ^F_ n )/_ n _ 2 

n=0 



We have 



F 2 (x,y) = -4, = 8(w(x) - 

16 Z" 2 

F (x,y) = -8(w(x)-w(y)) 2 , F_i(x, y) = — (u(x) - uj{y)f + A J J 2 



The local terms are 



(100) 



Lh 



o Jo 



dx I dy(X(x)Y{y) - X(y)Y(x)) { -S'"(x - y) + p 2 5'(x - y)+ 



-{uj{x) - uj{y))5"{x -y)- 2(oj(x) - oj{y))H'{x - y) - -F^x, y)5(x - y) 

The last two terms obviously vanish and what remains is 

L b = f dx{-(X"'Y - XY'") - A((X'Y - XY')(p 2 + 2a;')} 
Jo 



Lemma 7 The non local term eg. \W(J\) is identically zero. 



Proof . The non local term reads 

-l r i 



NL h 



-2 1 dx dy{X{x)Y{y) - X(y)Y(x))e(x - y) 
'0 Jo 

oo 



n=0 



£ (F^ n (x,y) + plF_ n (x,y))—(x-yT 



The first sum is just the coefficient of (iZ) in the formal expansion of dw (x, Z)dw* (y, Z) while 
the second sum is the coefficient of (iZ) . Hence we have 

NL b = 2 I"* dx I dy(X(x)Y(y) - X{y)Y{x))e{x - y) f ^-(Z - p^Z' 1 ) dw 2 (x, Z)dw* 2 (y, Z) 

J J *» Coo 

The above expression is zero in the following sense. Let us write 



^{Z-p 2 Z-^) dw 2 (x,Z)dw* 2 (y,Z) = £ (?/ ~ 



(7/ -> 



^-{Z-plZ' 1 ) dw 2 (x,Z)d l+1 w* 2 (x,Z) 



We will show that all the integrals around in the right hand side are identically zero. 
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Since the function w(x,Z) satisfies the Schroedinger equation, its square w 2 (x,Z) satisfies a third 
order differential equation 

Vw 2 {x, Z) = -4Z 2 dw 2 {x, Z), V = d 3 + 8to'd + V 

Let us introduce a pseudo differential operator $ such that 

w 2 (x,Z) = <$>e 2iZx 

Then 

V = d^d 2 ^ 1 

Since T> is anti self-adjoint, we also have 

V = -V* = $*~ 1 d 2 $*d 
it follows that w* 2 (x, Z) which is solution of 

{-V*)w* 2 (x,Z) = -4Z 2 dw* 2 (x,Z) 

can be written as 

w* 2 (x,Z) = d- 1 ^- 1 de~ 2iZx 

Hence 

dw 2 = d$e 2iZx , d i+1 w* 2 = d l §*~ l de- 2iZx 
Finally, we have to compute 

f —(Z -p^Z- 1 ) (dw 2 (x,Z))d i+1 w* 2 (x,Z) = [ —(Z-plZ-^idQe^^^de- 21 ^ 

JC^ 2&7T J Coo 2lTT 

= —[ — (d^e 2iZx )d i <f*- 1 (d 2 + Ap 2 )e- 2iZx 
2i lr 2m 

We recall the formula (see e.g. |15j) 

/ {De iZx ){Fe~ iZx ) = Res d (DF*) 

Jc x 2zvr 

where Res^ is Adler's residue So our expression is equal to 

Res d (d<f>{d 2 + 4pg)$~ 1 d i ) = Res d ({V + 4pgd)d*) = 
because (P + <ip 2 ] d)d % is a differential operator. ■ 

Consider next the term proportional to a(Z) in eq. (|98|) . 

4 f dx f dy[x(x)Y{y) - X(y)Y(x)) f ^-a(Z) (l - &\ d x w 2 (x, Z)d y (w(y, Z) w*(y, Z)) 

J J J Coo \ / 
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Lemma 8 Let us define 



One has 



dZ 



J P (x) = I ^a(Z)(iZfe 2lZx , p = -1, -2, • 

' Coo 



J_i(x) = ~6(x), J- p (x)=2P- 



rV-1 



(p-2) 



■(e(x)-l), p>2 



Proof . One has 

d x Jp(x) = 2J p+1 (x) 

The calculation of J_i(x) is easy. Next, we need boundary conditions to determine the other J p 
by integration 

- 2 r dZ e- lZ 1 



£(iO"J-„(0) = fjjf 
n=2 17 



2t 



2i7r sin Z Z — Q 

e~< x - 1 

+ / 7 

sm C C — 



2/ 



sin£ 



+ cotC 



2 



It follows that all the non local terms containing J p (x) for p < — 2 vanish when < a; < 1. The 
a(Z) term is 

L a = A I dx J dy(x(x)Y(y) - X(y)Y(x)^A 1 {x)C^ 2 (y)J^ 1 (x) 



or 



L a = 4[ dxX(x)5{x) [ dyY{y)J'(y)-A f dxY{x)5(x) [ dyX(y)io"(y) 
Jo Jo Jo Jo 

Finally, the term in c(Z) is just equal to the a(Z) one and double it. 

Putting everything together, we arrive at ea. ()93|) . In the course of this proof, we have shown the 
identities 

[ ^-(Z-pZZ- 1 ) (dw 2 (x,Z))d i+1 w* 2 (x,Z) = 0, Vi>0 

,J Woo 

These are quartic identities on the coefficients of w(x,Z). Of course, we also have the quadratic 
relations of Hirota and Sato that were interpreted by Sato as Plucker relations defining the infinite 
Grassmannian, allowing to give a precise definition of r-functions. The above relations are quartic 
relations on r-functions analogous to the quartic relations on Riemann Theta functions. 



12 Poisson bracket {Lo,u(y)}. 



In the previous section, we have obtained the Poisson bracket for the generators T n = L ri — Lq. We 
now have to reintroduce Lq and check that it has the correct Poisson brackets. The candidate for 
L was given in ea. (|66j) . Let us recall it: 

Lq=Pq + 
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where 
and 

, (Z) = ! - ^ _ ! _ ^ Sftf 

^Z 2 -m 2 7r 2 ^ Z 2 - m 2 7r 2 

m m 

where we have used ea, (|77|) . evaluated at x = 0, to express r\ m . 

Proposition 20 We have the following Poisson bracket: 

-1 r l 



L ,y o dyY(y)u(y) j = -4 jf d y y(y) [(L - pl)8'{y) + T'(y)j (101) 
This shows that {Lq, •} acts on u(y) as d y , as it should be. 

Again, the proof is long and we will split it into several lemmas. We need to compute {<f> m (x), T(y)}. 
and {<&' m (x), T(y)} for x = and x = 1. 

Multiplying ea. (|94|) by E n (y) and remembering that (i^G -1 = (F\ and F k {x) = — ^ F k {x), we 
get 

{$ m (a0, (£(y)*(jO>} = (102) 
= -2<t> m (x) £ Sm / fc 2 Cfc2 L~ 2 P _° Z 2 2) 2 , fe(x)F fe (x)F fe (y)G fc (y) - Fi(y)F>( X )G k (x) 

+2<(x) £ Sln / fc ?t~/° Z t \, \n^)h{y)G k (y) - Fl(y)F k (x)G k {x) 



As before, this can be expressed in terms of S7i and J7,2- We find 



{<M*),<£(l/)*(l/)>} = - ^ Res ±Zk -i-^ ^($ m (x)^0 2 - <(x)0 2 ) 



This formula is the starting point to begin the computation of {f](po), {E(y)&(y))}. Setting x = 
in ea. (|102j) . the term proportional to & m (x) vanishes because F k (0) = 0. We are left with (x = 0, 
but we keep it for a while) 



{$ m (x),(E(y)*(y))} = -$ m (x)^Res ±z ^ ffj^ (103) 

A; ^ 

Multiplying by , a m a 2^ an d summing over m, in the right hand side appears the sum 

(p —m it j 

^ £4(o)<Mo) 

VZ 2 7T 2 777 2 \ r f -TT 2 m 2 

Tr2 /vtq 2 ^ 



{ (pg - TO 2 7T 2 )(Z 2 - 7r 2 m 2 ) ^ (pfi - vr 2 m 2 )(Z 2 - vr 2 m 2 ) 

3^2 2^ 



Z| — j»q \Pq — m 2 7r 2 Z| — 7r 2 m 2 
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where we used that J2m z' 1 -™ 12 2 = ^' ^* ^ e factor 1/(Z% — Pq) cancels with the factor (1 

£ TV Ttl 

PqZ^ 2 ) in ea. p03[) and we are left with 

— tl ( Tin 1 \ RfiQ I ry 

: V 2 (Z)Z 2 



Finally 



V(Po),J q dyY{y)T(y)^=2r ] {pv) ^ dyY{y) jf ^-L^Af^) (104) 

where we used the asymptotics expressions for w and w* inside Q\ and Q2, hence removing the 
Tj 2 {Z) factor. The last step consists in evaluating the integral around C^. The reult is: 



Lemma 9 ^ 

»KPo),jf d 2/ y(y)T(y)| = Mpo) dyY(y)5'(y) (105) 

Proof . Let us consider the integral over in eq. (|104|) . The term containing b(Z) can be written 



MS 



Lb = L x {( g ( x ^ z )^*'( x ^ z )-^'( x ^ z )^*( x ^ z ))Mw(y^ z )^*(y^ z )) + 

Z)w(x, Z)d y w* 2 (y, Z) - w*(x, Z)w*'(x, Z)d y w 2 (y, Z)\ 



w(x, 



On the first line, we recognize the wronskian of w(x, Z) and w*(x, Z), which is just equal to —2iZ. 
Since d y (w(y, Z)w*(y, Z)) = Z~ 2 S' 2 (y) + • • • this term vanish by ea. (|99|) . Hence 



dZ b(Z) 



' x w 2 (x, Z)d y w* 2 (y, Z) - 8 x w* 2 (x, Z)d y w 2 (y, Z)} (106) 



2 Jr 2m Z 2 

OO 

= - F-.p(x,y)I- p -2(x - y) 

p=-2 

= 5'{y-x) + \e{x-y) jf ^I^ 2 (x, Z)^* 2 (y, Z) (107) 
The e(x — y) term is zero because 

f ^^ +1 w 2 (x,Z)dw* 2 (x,Z) = f ^Lhd^ x ^ Zx )^- x de- 2iZx ) 
Jr 2m Z Jr 2m Z 

,J Woo " 00 

Jr 2m 
= Res a (d m $ • = 



Next, we consider the a(Z) term. It reads 

dZ a(Z 



L a = 9 y I ^^^-(w(x,Z)w*(y,Z) - w*{x,Z)w(y,Z))w'(x,Z)w(y,Z) 
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dZ " {/) l/l d x w 2 (x,Z)d y (w(y,Z)w*(y,Z)) 



r 2in Z 2 I 2 

*s CO 



-(w*(x, Z)w'(x, Z) + w* (x, Z)w(x, Z))d y w 2 (y, Z) 
--(w*(x, Z)w'(x, Z) — w* (x, Z)w(x, Z))d y w 2 (y, Z) 



Again, the last term is the wronskian and so 

L a = [ ^^(-iZd y (w 2 (y,Z) + ^d x w 2 (x,Z)d y (w(y,Z)w*(y,Z)) 

" Coo V 



-\d x {w*(x, Z)w(x, Z))d y (w 2 (y, Z)) 



The first term is 



f ^-a(Z)^Z)-^ 1 A^ p ( y y iZy = E A^J^y) 

-'Coo p=—l 



i JUL (?v)p~ 1 

= -S'(y) - u>{y)5{y) + (e(y) - 1) £ A„ p (y)^— 

The last sum is zero because < y < 1. The 5{y) term vanishes because to(0) = 0. The second 
term is 

\ I ^^ A -P^)C- g (ym)-^e 2iZx = -he(x) - 1) (A 1 (x)C_ 2 (y)(2x) + • • • ) (108) 

J Coo 

This vanishes when x = 0. The third term is 

~\ [ ^^A^ p (y)C^(x)(zZ)-^e 2iZ y = \{e{y) - 1) (il 1 (y)CL. 2 (x)(2y) + • • •) 

J Coo 

and this vanishes when < y < 1. 

Finally, it is easy to see that the c(Z) term is equal to the a(Z) one. Putting everything together, 
we get ea. (|105f) . ■ 



The last result we need is the following: 
Lemma 10 

lm Po~™ -/o J A) 

Proof . Taking the derivative with respect to x of ea. ()lU2j) and remembering that F k (0) = Fj~(l) = 0, 
the remaining terms are (there is a cancellation in the & m (x) term) 

W m (x),(E(yMy)) } = -» B (^^ l ( !:ffi x 



5, 



^(xjFfcCxjFfc^Gfcd/) - F 2 (y)F' k (x)G k (x 
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where it is understood that x = or x = 1. By exactly the same argument as before 

^ Ef m (x)# {x) f 1 dyY{y)T{y)) = 2 (po) f 1 dyY (y)d y d x [ ^—i—flAO, 
V Po - n m Jo Jo Jc x 2ivr Z*rj(Z) 

Hence, we just have to take the derivative with respect to x of the previous result, before setting 
x = or x = 1. At x = we get 

27?(po) / ^y(y)[-<5"(y) + 4u/'( 2/ ) 
■/ o L 

The S"(y) term comes from eq. (|l()7jl while the second term comes from eq, (|108j) doubled by the c{Z) 
contribution. At x = 1 only the periodic b~"{y) remains. Taking the difference we obtain ea. ()lU9() . ■ 



Putting everything together, we arrive at ea. (|lUl|) . and this finishes our proof that u(x) does satisfy 
the Virasoro Poisson bracket. 



13 Conclusion. 



We have succeeded to take the continuum limit in the formulae expressing the dynamical variables 
of the Volterra model in terms of the separated variables. This yields exactly solvable potentials 
and formulae for the Virasoro generators of a rather unusual type. Still, we were able to check 
that they have the correct Poisson brackets. Of course the most interesting thing now is to try to 
quantize this approach. As a first step, a semiclassical analysis along the lines of |2H] should be 
very enlightening. The full quantum theory however may reserve some surprise. The bracket eq. (JSJ) 
being in fact an ordinary quadratic bracket, it is natural to quantize it with Weyl type commutation 
relations. This opens up the possibility of phenomena as those advocated in |§1 I27| . 
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